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Abstract 

The recently formulated theory of horizontal visibility graphs transforms time series into graphs 
and allows the possibility of studying dynamical systems through the characterization of their 
associated networks. This method leads to a natural graph-theoretical description of nonlinear 
systems with qualities in the spirit of symbolic dynamics. We support our claim via the case study 
of the period-doubling and band-splitting attractor cascades that characterize unimodal maps. 
We provide a universal analytical description of this classic scenario in terms of the horizontal 
visibility graphs associated with the dynamics within the attractors, that we call Feigenbaum 
graphs, independent of map nonlinearity or other particulars. We derive exact results for their 
degree distribution and related quantities, recast them in the context of the renormalization group 
and find that its fixed points coincide with those of network entropy optimization. Furthermore, 
we show that the network entropy mimics the Lyapunov exponent of the map independently of its 
sign, hinting at a Pesin-like relation equally valid out of chaos. 
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Introduction 

We expose a remarkable relationship between nonlinear dynamical systems and complex 
networks by means of the horizontal visibility (HV) algorithm [l-3] that transforms time 
series into graphs. In low- dimensional dissipative systems chaotic motion develops out of 
regular motion in a small number of ways or routes, and amongst which the period-doubling 
bifurcation cascade or Feigenbaum scenario is perhaps the better known and most famous 
mechanism This route to chaos appears an infinite number of times amongst the 

family of attractors spawned by unimodal maps within the so-called periodic windows that 
interrupt stretches of chaotic attractors. In the opposite direction, a route out of chaos 
accompanies each period-doubling cascade by a chaotic band-splitting cascade, and their 



shared bifurcation accumulation points 
known to possess universal properties 



brm transitions between order and chaos that are 
4|-[6). Low- dimensional maps have been extensively 



studied from a purely theoretical perspective, but systems with many degrees of freedom 
used to study diverse problems in physics, biology, chemistry, engineering, and social 
science, are known to display low- dimensional dynamics 3]. 



The horizontal visibility (HV) algorithm converts the information stored in a time series 



into a network, setting the nature o: 



requires complex network tools 



the dynamical system into a different context that 



-Qj| to extract its properties. This approach belongs to 



an emerging corpus of methods that map series to networks (see for instance 2|, Il3l-ll7j or 
a recent review |18|). Relevant information can be obtained through the family of visibility 
methods, including the characterization of fractal behavior [l9[ or the discrimination 
between random and chaotic series [l. 
fields, from geophysics |2lJ, to finance 
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and it finds increasing applications in separate 
or physiology J^. Here we offer a distinct view 
of the Feigenbaum scenario through the specific HV formalism, and provide a complete set 
of graphs, which we call Feigenbaum graphs, that encode the dynamics of all stationary 
trajectories of unimodal maps. We first characterize their topology via the order-of-visit 
and self-affinity properties of the maps. Additionally, a matching renormalization group 
(RG) procedure leads, via its flows, to or from network fixed-points to a comprehensive 
view of the entire family of attractors. Furthermore, the optimization of the entropy 
obtained from the degree distribution coincides with the RG fixed points and reproduces 
the essential features of the map's Lyapunov exponent independently of its sign. A general 
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observation is that the HV algorithm extracts only universal elements of the dynamics, 
free of the peculiarities of the individual unimodal map, but also of universality classes 
characterized by the degree of nonlinearity. Therefore all the results presented in this work, 
while referring to the specific Logistic map for illustrative reasons apply to any unimodal 
map. 



Model: Feigenbaum graphs 

The HV graph [1] associated with a given time series {xi}i=i,...,jv of iV real data is 
constructed as follows: First, a node % is assigned to each datum Xj, and then two nodes 

' a,ld 3 m co,mected if the corre8poildmg data fulfi " the > x - ** a " " 

such that i < n < j. Let us now focus on the Logistic map [4J defined by the quadratic 
difference equation x t +\ = f(xt) = ^x t (l — x t ) where x t G [0, 1] and the control parameter 
G [0,4]. According to the HV algorithm, a time series generated by the Logistic map for 
a specific value of \l (after an initial transient of approach to the attractor) is converted 
into a Feigenbaum graph (see figured]). Notice that this is a well-defined subclass of HV 
graphs where consecutive nodes of degree k = 2, that is, consecutive data with the same 
value, do not appear, what is actually the case for series extracted from maps (besides the 
trivial case of a constant series). While for a period T there are in principle several possible 
periodic orbits, and therefore the set of associated Feigenbaum graphs is degenerate, it 
can be proved that the mean degree k(T) and normalized mean distance d(T) of all these 
Feigenbaum graphs fulfill k(T) = 4(1 — ^) and d{T) = respectively, yielding a linear 
relation d{k) = (4 — k)/6 that is corroborated in the inset of figure [Q Observe that aperiodic 
series (T — > oo) reach the upper bound mean degree k = 4. 



Results 

A deep-seated feature of the period-doubling cascade is that the order in which the positions 
of a periodic attractor are visited is universal 24j . the same for all unimodal maps. This 
ordering turns out to be a decisive property in the derivation of the structure of the Feigen- 
baum graphs. See figure [2] where we plot the graphs for a family of attractors of increasing 
period T = 2 n , that is, for increasing values of fi < fi^. This basic pattern also leads to the 
expression for their associated degree distributions, 



4 



P(n,k) = {lf\ fc = 2,4,6,...,2n, (1) 
P(n,k) = Q) n , fc = 2(n + l), 

and zero for k odd or > 2(n + 1). At the accumulation point /Xqq the period diverges 
(n oo) and the distribution is exponential for all even values of the degree, 



F(oo,A;)=f-l , fc = 2,4,6,..., (2) 

and zero for odd. Observe that these relations are independent of the order of the map's 
nonlinearity: the HV algorithm sifts out every detail of the dynamics except for the basic 
storyline. 

We turn next to the period-doubling bifurcation cascade of chaotic bands that takes place 
as fi decreases from /i = 4 towards fi^. For the largest value of the control parameter, at 
H — 4, the attractor is fully chaotic and occupies the entire interval [0, 1] (see figured]). This 
is the first chaotic band n = at its maximum amplitude. As fi decreases in value within 
/^oo < /i < 4 band-narrowing and successive band-splittings occur. In general, after 



n reverse bifurcations the phase space is partitioned in 2 n disconnected chaotic bands, which 



are self-affine copies of the first chaotic band 25j. The values of /x at which the bands split 
are called Misiurewicz points 24], and their location converges to the accumulation point 
fioo for n — > oo. Significantly, while in the chaotic zone orbits are aperiodic, for reasons 
of continuity they visit each of the 2 n chaotic bands in the same order as positions are 
visited in the attractors of period T = 2" [24]. In figure [3] we have plotted the Feigenbaum 
graphs generated through chaotic time series at different values of fi that correspond to 
an increasing number of reverse bifurcations. Since chaotic bands do not overlap, one can 
derive the following degree distribution for a Feigenbaum graph after n chaotic-band reverse 
bifurcations by using only the universal order of visits 

P ll (n,k)= (|) fc/2 , fc = 2,4,6,...,2n, 

P,(n,k>2(n + 1))= (i) n , (3) 

and zero for k = 3,5, 7, 2n+l. We note that this time the degree distribution retains some 
dependence on the specific value of ji, concretely, for those nodes with degree k > 2(n + 1), 



all of which belong to the top chaotic band (labelled with red links in figure |3]). The HV 
algorithm filters out chaotic motion within all bands except for that taking place in the top 
band whose contribution decreases as n — » oo and appears coarse-grained in the cumulative 
distribution P^(n,k > 2(n + 1)). As would be expected, at the accumulation point fi^ 
we recover the exponential degree distribution (equation [2]), i.e. lim n ^oo P^{n, k) = P(oo, k). 

Before proceeding to interpret these findings via the consideration of renormalization 
group (RG) arguments, we recall that the Feigenbaum tree shows a rich self-affine structure: 
for fi > /ioo periodic windows of initial period m undergo successive period-doubling 
bifurcations with new accumulation points /^(m) that appear interwoven with chaotic 
attractors. These cascades are self-affine copies of the fundamental one. The process of re- 
verse bifurcations also evidences this self-affine structure, such that each accumulation point 
is the limit of a chaotic-band reverse bifurcation cascade. Accordingly, we label G(m, n) 
the Feigenbaum graph associated with a periodic series of period T = m ■ 2 n , that is, a 
graph obtained from an attractor within window of initial period m after n period-doubling 
bifurcations. In the same fashion, G M (n, m) is associated with a chaotic attractor composed 
by m ■ 2 n bands (that is, after n chaotic band reverse bifurcations of m initial chaotic 
bands). Therefore, graphs depicted in figures [2] and [3] correspond to G(l,n) and G^(l,n) 
respectively and for the first accumulation point we have G(l,oo) = 6^(1,00) = G^. 
Similarly, in each accumulation point /ioo(^) we have G(m, 00) = G M (m, 00). 



In order to recast previous findings in the context of the renormalization group, let us 
define an RG operation TZ on a graph as the coarse-graining of every couple of adjacent 
nodes where one of them has degree k = 2 into a block node that inherits the links of 
the previous two nodes (see figure Hla). This is a real-space RG transformation on the 
Feigenbaum graph [261 ] 
renormalization schemes 



aissimi 



ar from recently suggested box-covering complex network 
27M29I]. This scheme turns out to be equivalent for /i < /i^ to 
the construction of an HV graph from the composed map instead of the original /, in 
correspondence to the original Feigenbaum renormalization procedure 0, [3^] . We first note 
that TZ{G(l,n)} = G(l,n — 1), thus, an iteration of this process yields an RG flow that 
converges to the (1st) trivial fixed point TZ^{G{l,n)} = G(1,0) = G Q = TZ{G }. This is 
the stable fixed point of the RG flow V/i < /x^. We note that there is only one relevant 
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variable in our RG scheme, represented by the reduced control parameter A/x = fi^ — /i, 
hence, to identify a nontrivial fixed point we set A/i = or equivalently n — > oo, where the 
structure of the Feigenbaum graph turns to be completely self-similar under TZ. Therefore 
we conclude that G(l, oo) = G^ is the nontrivial fixed point of the RG flow, TZ{G oc } = G^. 
In connection with this, let Pt(k) be the degree distribution of a generic Feigenbaum graph 
Gt in the period-doubling cascade after t iterations of TZ, and point out that the RG 
operation, TZ{G t } = G t +i, implies a recurrence relation (1 — P t {2))P t+ i{k) = P t {k + 2), 
whose fixed point coincides with the degree distribution found in equation [2j This confirms 
that the nontrivial fixed point of the flow is indeed G^. 

Next, under the same RG transformation, the self-affine structure of the family of 
attractors yields TZ{G fl (l,n)} = G^l, n — 1), generating a RG flow that converges to 
the Feigenbaum graph associated to the 1st chaotic band, TZ^{G^(l,n)} = G^(1,0). 
Repeated application of TZ breaks temporal correlations in the series, and the RG flow 
leads to a 2nd trivial fixed point TZ^lG^l, 0)} = G ran d = TZ{G ISkn d}, where G ran d is the 
HV graph generated by a purely uncorrelated random process. This graph has a universal 
degree distribution P(k) = (l/3)(2/3) fc ~ 2 , independent of the random process underlying 
probability density (see jl[ Q). 

Finally, let us consider the RG flow inside a given periodic window of initial period m. 
As the renormalization process addresses nodes with degree k = 2, the initial applications 
of TZ only change the core structure of the graph associated with the specific value m 
(see figure Hlb for an illustrative example). The RG flow will therefore converge to the 
1st trivial fixed point via the initial path TZ^{G(m, n)} = G(l,n), with p < m, whereas 
it converges to the 2nd trivial fixed point for G^(m,n) via TZ^ {G ^(m, n)} = G^(l, n). 
In the limit of n — > oo the RG flow proceeds towards the nontrivial fixed point via the 
path TZ^{G(m, oo)} = G(l,oo). Incidentally, extending the definition of the reduced 
control parameter to A/i(m) = /^(m) — /i, the family of accumulation points is found 
at A/x(m) = 0. A complete schematic representation of the RG flows can be seen in figure Hlc. 

Interestingly, and at odds with standard RG applications to (asymptotically) scale- 
invariant systems, we find that invariance at A/z = is associated in this instance to an 
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exponential (rather than power-law) function of the observables, concretely, that for the 
degree distribution. The reason is straightforward: 1Z is not a conformal transformation (i.e. 
a scale operation) as in the typical RG, but rather, a translation procedure. The associated 
invariant functions are therefore non homogeneous (with the property g(ax) = bg(x)), but 
exponential (with the property g(x + a) = cg(x)). 

Finally, we derive, via optimization of an entropic functional for the Feigenbaum 
graphs, all the RG flow directions and fixed points directly from the information contained 
in the degree distribution. Amongst the graph theoretical entropies that have been 
proposed we employ here the Shannon entropy of the degree distribution P(k), that is 
h = — ^2^=2 log-P(Zc). By making use of the Maximum Entropy formalism, it is easy to 
prove that the degree distribution P(k) that maximizes h is exactly P(k) = (l/3)(2/3) fc_2 , 
which corresponds to the distribution for the 2nd trivial fixed point of the RG flow Grand- 
Alternatively, with the incorporation of the additional constraint that allows only even 
values for the degree (the topological restriction for Feigenbaum graphs G(l,n)), entropy 
maximization yields a degree distribution that coincides with equation El which corresponds 
to the nontrivial fixed point of the RG flow G^. Lastly, the degree distribution that 
minimizes h trivially corresponds to Go, i-e. the 1st trivial fixed point of the RG flow. 
Remarkably, these results indicate that the fixed-point structure of the RG flow are obtained 
via optimization of the entropy for the entire family of networks, supporting a suggested 
connection between RG theory and the principle of Maximum Entropy j3ll |. 

The network entropy h can be calculated exactly for G(l,n) (/i < floo or T = 2 n ), yielding 
h(n) = log 4 • (1 — 2 _n ). Because increments of entropy are only due to the occurrence of 
bifurcations h increases with fi in a step-wise way, and reaches asymptotically the value 
h{po) = log 4 at the accumulation point /Xqq. For Feigenbaum graphs G M (1, n) (in the chaotic 
region), in general h cannot be derived exactly since the precise shape of P(k) is unknown 
(albeit the asymptotic shape is also exponential [20j). Yet, the main feature of h can be 
determined along the chaotic-band splitting process, as each reverse bifurcation generates 
two self-affine copies of each chaotic band. Accordingly, the decrease of entropy associated 
with this reverse bifurcation process can be described as h^(n) = log4 + h fl (0)/2 n , where 
the entropy h^in) after n reverse bifurcations can be described in terms of the entropy 
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associated with the first chaotic band fyu(O). In figured] we observe how the chaotic-band 
reverse bifurcation process takes place in the chaotic region from right to left, and therefore 
leads in this case to a decrease of entropy with an asymptotic value of log 4 for n — > oo at 
the accumulation point. These results suggest that the graph entropy behaves qualitatively 
as the map's Lyapunov exponent A, with the peculiarity of having a shift of log 4, as 
confirmed in figure [5j This unexpected qualitative agreement is reasonable in the chaotic 
region in view of the Pesin theorem j^], that relates the positive Lyapunov exponents of a 
map with its Kolmogorov-Sinai entropy (akin to a topological entropy) that for unimodal 
maps reads hxs = A, VA > 0, since h can be understood as a proxy for hxs- Unexpectedly, 
this qualitative agreement seems also valid in the periodic windows (A < 0), since the 
graph entropy is positive and varies with the value of the associated (negative) Lyapunov 
exponent even though hxs — 0, hinting at a Pesin-like relation valid also out of chaos which 
deserves further investigation. The agreement between both quantities lead us to conclude 
that the Feigenbaum graphs capture not only the period-doubling route to chaos in a 
universal way, but also inherits the main feature of chaos, i.e. sensitivity to initial conditions. 



Discussion 

In summary, we have shown that the horizontal visibility theory combines power with 
straightforwardness as a tool for the analytical study of nonlinear dynamics. As an illus- 
tration we have established how the families of periodic and chaotic attractor bifurcation 
cascades of unimodal maps transform into families of networks with scale-invariant limiting 
forms, whose characterization can be deduced from two basic and universal properties of 
unimodal maps: ordering of consecutive positions in the attractors and self-affinity Further, 
we have demonstrated that these networks and their associated degree distributions comply 
with renormalization group and maximum entropy principles, filtering out irrelevant 
variables and finding fixed-point networks which are independent of the map's nonlinearity. 
The entire Feigenbaum scenario is therefore fully described. The potential of the theory for 
revealing new information is indicated by the ability of the network entropy to emulate the 
Lyapunov exponent for both periodic and chaotic attractors. Extensions of this approach 
to other complex behavior, such as dynamical complexity associated to vanishing Lyapunov 
exponents, intermittency, quasiperiodic routes to chaos, etc., are still open questions. 
Finally, observe that in symbolic dynamics 32( one usually defines a phase-space partition 
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(Markov partition) in order to create a symbolic representation of the dynamics. This 
partition tiles phase space in a non-overlapping manner: every value of the series has a 
univocally associated symbol. While a Feigenbaum graph also symbolizes the series data 
(incidentally, without the need of defining an ad hoc partition), each series datum is not 
associated univocally to a symbol (the degree of the node): this symbol is a function, in 
principle, of the complete series, and incorporates global information. Furthermore, note 
that besides the symbolization that converts a time series into a series of node degrees, the 
Feigenbaum graphs also store the connectivity pattern amongst nodes -i.e. the topological 
structure of the graph. On this respect, the possible connections of HV theory with 
kneading theory |33| and symbolic dynamics [32| are of special interest. 
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Figure Legends 
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FIG. 1: Feigenbaum graphs from the Logistic map x t +i = f(xt) = [ix t (l — xt). The main 
figure portrays the family of attractors of the Logistic map and indicates a transition from periodic 
to chaotic behavior at = 3.569946... through period-doubling bifurcations. For \i > fi^ the 
figure shows merging of chaotic-band attractors where aperiodic behavior appears interrupted by 
windows that, when entered from their left-hand side, display periodic motion of period T = m ■ 2° 
with m > 1 (for /i < fj,^, m = 1) that subsequently develops into m period-doubling cascades 
with new accumulation points /^(m). Each accumulation point fi^fn) is in turn the limit of a 
chaotic-band reverse bifurcation cascade with m initial chaotic bands, reminiscent of the self-affine 
structure of the entire diagram. All unimodal maps exhibit a period-doubling route to chaos with 
universal asymptotic scaling ratios^between successive bifurcations that depend only on the order 
of the nonlinearity of the map 



30l | . the Logistic map belongs to the quadratic case. Adjoining 
the main figure, we show time series and their associated Feigenbaum graphs according to the 
HV mapping criterion for several values of /i where the map evidences both regular and chaotic 
behavior (see the text). Inset: Numerical values of the mean normalized distance d as a function 
of mean degree k of the Feigenbaum graphs for^i < /i < 4 (associated to time series of 1500 data 
after a transient and a step 5fi = 0.05), in good agreement with the theoretical linear relation (see 




FIG. 2: Periodic Feigenbaum graphs for \i < /j,^. The sequence of graphs associated to 
periodic attractors with increasing period T = 2 n undergoing a period-doubling cascade. The 
pattern that occurs for increasing values of the period is related to the universal ordering with 
which an orbit visits the points of the attractor. Observe that the hierarchical self-similarity of 
these graphs requires that the graph for n — 1 is a subgraph of that for n. 
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u =3.5748... 

M-8 

FIG. 3: Aperiodic Feigenbaum graphs for fi > fj,^. A sequence of graphs associated with 
chaotic series after n chaotic-band reverse bifurcations, starting at fj, = 4 for n = 0, when the 
attractor extends along a single band and the degree distribution does not present any regularity 
(red links). For n > the phase space is partitioned in 2 n disconnected chaotic bands and the 
n-th self-afnne image of [i = 4 is the n-th Misiurewicz point /i 2 n - 1 -2 n - ln an cases, the orbit visits 
each chaotic band in the same order as in the periodic region [i < fi^. This order of visits induces 
an ordered structure in the graphs (black links) analogous to that found for the period-doubling 
cascade. 
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FIG. 4: Renormalization process and network RG flow structure, (a) Illustration of the 
renormalization process 1Z: a node with degree k = 2 is coarse-grained with one of its neighbors 
(indistinctively) into a block node that inherits the links of both nodes. This process coarse-grains 
every node with degree k = 2 present at each renormalization step, (b) Example of an iterated 
renormalization process in a sample Feigenbaum graph at a periodic window with initial period 
m = 9 after n = 2 period-doubling bifurcations (an orbit of period T = m ■ 2 n = 36). (c) RG flow 
diagram, where m identifies the periodic window that is initiated with period m and n designates 
the order of the bifurcation, n = n + 1 for period-doubling bifurcations and n = — (n + 1) for 
reverse bifurcations. A / u(m) = ji^m) — [i denotes the reduced control parameter of the map, and 
V-ooijn) is the location of the accumulation point of the bifurcation cascades within that window. 
Feigenbaum graphs associated with periodic series (A/i(m) > 0, n > 0) converge to (5(1,0) = Go 
under the RG, whereas those associated with aperiodic ones (A/i(m) < 0, n < 0) converge to 
Grand- The accumulation point = /ioo(l) corresponds to the unstable (nontrivial) fixed point 
G(l,oo) = Gqo of the RG flow, which is nonetheless approached through the critical manifold of 
graphs G(m, oo) at the accumulation points /^(m). In summary, the nontrivial fixed point of the 
RG flow is only reached via the family of the accumulation points, otherwise the flow converges to 
trivial fixed points for periodic or chaotic regions. 
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FIG. 5: Horizontal visibility network entropy h and Lyapunov exponent A for the 
Logistic map. We plot the numerical values of h and A for 3.5 < fi < 4 (the numerical step 
is 5/j, = 5 • 10~ 4 and in each case the processed time series have a size of 2 12 data). The inset 
reproduces the same data but with a rescaled entropy h — log(4). The surprisingly good match 
between both quantities is reminiscent of the Pesin identity (see text). Unexpectedly, the Lyapunov 
exponent within the periodic windows (A < inside the chaotic region) is also well captured by h. 
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